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ABSTRACT 



We analyze new and published optical photometric data of FU Orionis, an 
eruptive pre-main-sequence star. The outburst consists of a 5.5 mag rise at B 
with an e-folding timescale of ~ 50 days. The rates of decline at B and V are 
identical, 0.015±0.001 mag yr _1 . Random fluctuations superimposed on this 
decline have an amplitude of 0.035 ± 0.005 mag at V and occur on timescales 
of 1 day or less. Correlations between V and the color indices U-B, B-V, and 
V-R indicate that the variable source has the optical colors of a GO supergiant. 
We associate this behavior with small amplitude flickering of the inner accretion 
disk. 



Subject headings: accretion disks - stars pre-main-sequence - stars: formation 
stars: individual (FU Orionis) 
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1. INTRODUCTION 



FU Orionis objects - sometimes known as FUors - are eruptive pre-main-sequence 
stars located in active star forming regions ( [Herbig 1966| , 1977; |Hartmann fc Kenyon 199"E 



Kcnyon 1999 ). Roughly half of the 11 commonly accepted FUors have been observed to rise 
3-5 mag in optical or near-IR brightness on timescales of 1-10 yr. Other FUors have been 
identified based on properties similar to eruptive FUors, including (i) absorption features 



of F-G supergiants on optical spectra and K-M giants on near-IR spectra flHerbig 1977 



Mould et a/.T978[ |Carr et al. 19871 ; [Stocke et aZ.~1988l [Staude fc Neckel 1992|) ; (ii) large 



excesses of radiation over normal F-G stars at ultraviolet, infrared, sub millimeter, and 
centimeter wavelengths flWcintraub et al. 1989| , 1991; Kenyon fc Hartmann 199l| ; [Rodriguez 



et al. 1990| ; [Rodriguez fc Hartmann 1992| ); (iii) distinctive reflection nebulae ( [Goodrich 



|1987[ ); and (iv) clear association with optical jets, HH objects, and molecular outflows 



( [Uoipurth 1991 ; l^vans et al. 1994| ). 



FUor eruptions are often accepted as accretion events in the disk surrounding a 
low-mass pre-main-sequence star (|Hartmann fc Kenyon 19"85| ; |Lin fc Papaloizou 1985 



Hartmann fc Kenyon 1996| ; for alternative interpretations, see Herbig & Petrov 1992; 



Petrov et al. 1998). In this picture, the accretion rate through the disk increases by 2-3 
orders of magnitude to ~ 10~ 4 M Q yr -1 . In addition to providing energy for the luminosity 
increases of FUors, this model naturally explains the broad spectral energy distributions, 
the variation of rotational velocity and spectral type with wavelength, and color changes 
during the optical decline of V1057 Cyg, among other observed properties. 

Despite the success of the disk hypothesis, one observable characteristic of disk 
accreting systems - flickering - has not been observed in any known FUor. In most systems 
with luminous accretion disks, flickering is observed as a series of random brightness 
fluctuations with amplitudes of 0.01-1.0 mag that recur on dynamical timescales ( |Robinso~n 



1976; Bruch 1992| ). Often accepted as a 'signature' of disk accretion, flickering is believed 



to be a dynamical variation of the energy output from the diskF]. Thus, it provides some 
measure of the fluctuations in the physical structure of the disk and might someday serve 
as a diagnostic of physical properties within the disk ([Bruch fc Duschl 1993| ; [Bruch 19~94[ ; 



-'in the cataclysmic variables (short period binary systems with an accretion disk surrounding a white 
dwarf), random flickering occurs on timescales of seconds to minutes. Dwarf nova oscillations and quasi- 
periodic oscillations are semi-coherent periodic variations observed on similar timescales. Flickering is also 



occasionally associated with material in the 'bright spot' at the outer edge of the disk (see Warner 1995). 
It is unclear whether or not these classes have distinct analogs among accreting pre-main sequence stars. In 
this paper, we use flickering to distinguish rapid variations of light from the inner disk from variations of the 
bright spot. 
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Warner 1995|) . 

In this paper, we search for evidence of flickering in the historical light curve of FU 
Ori. In addition to the outburst, we find good evidence for small-amplitude brightness 
fluctuations on a timescale of 1 day or less. Color changes correlated with the brightness 
changes indicate a variable source with the optical colors of a GO supergiant. The amplitude, 
color temperature, and timescale of the variations have much in common with the flickering 
observed in short period interacting binary systems. After ruling out several possible 
alternatives, we conclude that flickering is the most likely interpretation for short-term 
variability in FU Ori. The most plausible location for the flickering source is the inner edge 
of the disk, where temperatures lie between the stellar temperature and the maximum disk 
temperature of ~ 7000 K. 

We describe the observations in §2, analyze the light curve in §3, and conclude with a 
brief discussion and summary in §4. 



2. OBSERVATIONS 

We acquired UBV photometry of FU Ori with the 60-cm Zeiss reflector at the Crimean 
Laboratory of the Sternberg State Astronomical Institute. The observations usually were 
made through a 13" aperture; a 27" aperture was used on nights of poor seeing. These data 
were reduced using BD+8°1051 as the comparison star and other nearby stars as controls 

Kolotilov fc Petrov 1985|). Table 1 lists the results. The uncertainty in the calibration 



[see 

is ±0.01-0.02 mag for V and B-V, and ±0.02-0.04 mag for U-B. 

We supplement the UBV data with additional photoelectric photometry from the 
Maidanak High Altitude Observatory. Ibragimov (1997) describes UBVR data acquired 
during 1981-1994 as part of the ROTOR project. The data have been reduced to the 
standard UBVR system with typical errors of ±0.015 mag for V and V-R, ±0.02 mag for 
B-V, and ±0.04-0.08 mag for U-B. 

We also consider visual observations of FU Ori compiled by the American Association 
of Variable Star Observers. The error of a typical estimate is 0.1-0.2 mag using standard 
stars in the field of FU Ori calibrated from photoelectric observations. We simplify a 
comparison with photoelectric data by computing twenty day means of the over 7000 
AAVSO observations. With ~ 15 observations per twenty day interval, this procedure 
reduces the typical error of a twenty day mean to ~ 0.03 mag, only slightly larger than the 
quoted error of the photoelectric data. 
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Finally, we obtained low resolution optical spectra of FU Ori during 1995-1998 with 
FAST, a high throughput, slit spectrograph mounted at the Fred L. Whipple Observatory 
1.5-m telescope on Mount Hopkins, Arizona QFabricant et at. 1998| ). We used a 300 g mm -1 
grating blazed at 4750 A, a 3" slit, and recorded the spectra on a thinned Loral 512 x 2688 
CCD. These spectra cover 3800-7500 A at a resolution of ~ 6 A. On photometric nights, 
we acquired standard star observations to reduce the FU Ori data to the Hayes & Latham 
(1975) flux scale using NOAO IRAFf]. These calibrations have an accuracy of ±0.05 mag. 
This uncertainty is comparable to the probable error in spectrophotometric data acquired 
with the Kitt Peak National Observatory Intensified Reticon Spectrograph reported in 
Kenyon et al. (1988). 

Table 2 lists indices for several absorption lines using O'Connell's (1973) definition, 
I\ = F\/F\, where F\ is the measured flux in a 20 A bandpass centered at wavelength A 
and F\ is the continuum flux interpolated from continuum bandpasses on either side of the 
absorption feature. Repeat measurements indicate an error of ± 0.03 mag for each index. 



3. LIGHT CURVE ANALYSIS 

The lower panel of Figure 1 shows historical B+photographic and optical light curves 
for FU Ori. Following a 5-6 mag rise at B, the system has declined by ~ 1 mag in nearly 
70 yr. The decline in visual light has closely followed the B light curve for the past 30 yr. 
In addition to a long-term wave-like variation, the brightness shows considerable scatter 
of roughly 0.1 mag at BVR and almost 0.2 mag in U on timescales of days to months. 
These fluctuations are much larger than the quoted photometric errors. Ibragimov (1997) 
commented on fluctuations in brightness and color indices throughout 1981-1994. Similar 
variations are visible in light curves of V1057 Cyg and V1515 Cyg ( [Ibragimov 1997] ). We 
plan analyses of these other FUors in future publications. 



3.1. Outburst Model 

To analyze brightness fluctuations in FU Ori, we begin by separating long-term changes 
with timescales of years from shorter timescale variations. We consider a simple model for 



2 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under contract to the National Science 
Foundation 
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the outburst 



m = m + m rise (t - t ) (1) 

where m is the average quiescent magnitude and t is the time of the outburst. The 
brightness during outburst is a function of time, t, 

m -1° t<t0 (2) 

nse \ 5m rise {l - e-^-^/— } + rh(t-t ) t > t {) 

where 5m r i se is the amplitude of the outburst, r r i se is the e- folding time of the rise to 
maximum, and m is the rate of decline from maximum. We derive parameters for this 
model using a downhill simplex method to minimize the function 

N 

x 2 = J2[(m-m t )/a t }\ (3) 

for N observations rrii having uncertainty 0{. We derive the model parameters mo, 5m r i se , 
rh, and r r j se for an adopted to and vary to separately to produce a minimum in \ 2 ■ This 
procedure works better than fitting all parameters at once due to the sparse nature of the 
light curve for t <^ t . To estimate errors of the model parameters, we compute residuals 
r,i of the light curve about the fit, construct new light curves by adding gaussian noise 
with amplitude r^ to the data, and extracting new model parameters. We adopt as best 
the median values for model parameters from 10,000 such trials; the quoted errors are the 
inter-quartile ranges of each model parameter. 

Table 3 summarizes results of model fits to the light curves in Figure 1. The sparse 
data before outburst yield a poor measure of the pre-outburst brightness and the start of 
the outburst. We estimate t = JD 2428497 ± 40 and m = 15.55 ± 0.40 from the B-band 
data; these errors set the uncertainties in the other model parameters. The uncertainty in 
the decline rate at B is small; the visual data yield a nearly identical rate of decline despite 
lack of data near maximum. 

The spectroscopic data show little evidence for any change in the mean optical spectral 
type during the recent 0.3 mag decline in the mean V brightness. The Mg I index, which 
tracks the spectral type rather well, increased by at most 0.02 mag during these years. 
The H/9 and Na I indices, which should measure the wind of FU Ori, have also remained 
constant. We show below that the mean colors, U-B, B-V and V-R, of FU Ori have 
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changed by <^ 0.02 mag in 15 years. The constant colors and spectral indices indicate that 
the optical spectral type has remained constant to within ± 1 subclass since 1985. 

The upper panel of Figure 1 shows residual light curves about the best-fit model. The 
rms dispersion about the light curves is 0.06 mag for photoelectric data, 0.11 mag for visual 
data, and 0.24 mag for photographic data. All are large compared to typical uncertainties. 

Some of the scatter in the residual light curves is due to a long-term wave with an 
amplitude of several tenths of a magnitude. The residual B light curve has three prominent 
crests at ~ JD 2428500, JD 2432000, and JD 2436000, with a possible small crest at ~ JD 
2433500. This behavior appears to vanish in the B light curve at later times. However, the 
residual visual light curve then shows a similar wave-like oscillation with crests at ~ JD 
2439000 and JD 2442500. Another crest at ~ JD 2447000 is present in the photoelectric V 
light curve shown in Figure 2 and described below. These peaks are roughly in phase with 
the residual B light curve for periods of ~ 4000 days. 

There are considerable photometric variations in addition to the wave. These occur 
on short timescales, tens of days or less, and have amplitudes that seem uncorrelated with 
the overall system brightness. We will analyze these in the next section and then consider 
possible origins for the short-term and long-term variations. 

Despite the large residuals, the model fits leave negligible trends in the data. The 
median Spearman rank correlation coefficient between the V brightness and time for the 
10,000 Monte Carlo trials is 0.87 for the B + photographic data and 0.84 for the visual 
data. The median slope of a linear fit to the residual light curves is ^ 10~ 4 mag yr _1 for 
the B + photographic data and ^ 10~ 6 mag yr _1 for the visual data. We conclude that the 
simple model provides a reasonable fit to the long-term light curve and now consider the 
nature of the shorter timescale fluctuations. 



3.2. Fluctuation Timescale 

To test whether or not the short timescale photometric variations in FU Ori are 
flickering, we need to verify that (i) they are real changes in the brightness and color of 
FU Ori, (ii) they occur on short timescales but are not periodic, and (iii) they can be 
plausibly associated with radiation from the disk. We establish the reality of the variations 
in two parts. We first demonstrate that the variations occur on short timescales with 
amplitudes larger than the photometric errors. We show later that color changes correlate 
with brightness changes. 
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We consider a Monte Carlo model for the photoelectric data described in §2. The lower 
panel of Figure 2 shows V-band data acquired during the last 15 yr. Both the wave-like 
oscillation and the large scatter about this oscillation are apparent. The solid line in this 
Figure is the seasonal mean light curve, the average brightness for each year of observation. 
The top panel in the Figure is the difference between the actual data and the seasonal mean. 
This residual light curve has a large amplitude but no linear trend or wave-like feature. 

Our analysis indicates a small periodic component in residual light curves for B and 
V data. Periodograms suggest a period of 17 ± 1 days in the B data, which has been 
reported previously ( [Kolotilov fc Petrov 1985|) . This periodic component has an amplitude 
of 0.009 ± 0.003 mag. The V data has a best period of 111.4 ± 1.6 days with an amplitude 
of 0.012 ± 0.003 mag. There is no indication of the 17 day period in the V light curve. 
The amplitudes of these 'periodicities' are comparable to the photometric errors but small 
compared to the amplitude of the fluctuations in the residual light curves (see Figures 1-2). 

To measure the amplitude of the non-periodic component in the residual light curve, 
we use a Monte Carlo model. We replace each observation Vj with a random brightness 
having amplitude a v and offset t>o, 

Vi = v + a v g { , (4) 

where gi is a normally distributed deviate with zero mean and unit variance ( Press et al. 



1992|) . Artificial light curves that provide a good match to the actual light curve should have 
the same amplitude and offset. We quantify a 'good match' by comparing the magnitude 
distributions of actual and artificial light curves using the Kolmogorov-Smirnov (K-S) test. 
We reject poor matches with a low probability of being drawn randomly from the same 
distribution as the data. The 'best' match maximizes the median K-S probability from 
10,000 trials. We establish error estimates for the best parameters and the proper scale for 
this measure by comparing two artificial light curves generated in each of 10,000 trials. 

This procedure yields best parameters of v = 0.0016±0.0009 mag and a v = 0.033±0.005 
mag for the residual V light curve. Artificial light curves with these parameters have a 
high probability, 68% or larger, of being drawn randomly from the same distribution as the 
actual data. The offset of the model light curve is consistent with zero. The amplitude of 
the model is roughly twice the quoted la error of 0.015 mag. 

The artificial data sets have periodic variations similar to the real data sets. Periods 
of 10-100 days are common in the artificial B and V light curves. Mean light curves folded 
on these periods have amplitudes, 0.01 ± 0.003 mag, similar to those quoted above for 
the real data. The 'best' period is different in each artificial data set, but the amplitude 



is nearly constant. This amplitude is small compared to the random component of the 
fluctuation. We suspect the periodic variations are due to the sampling of the light curve. 
Such 'periodicities' illustrate some of the dangers of period analysis. 

There are two explanations for the 0.03 mag variations in the residual V light curve. A 
simple test should distinguish real fluctuations from measurement error. Real fluctuations 
in V should be accompanied by correlated variations in the color indices, U-B, B-V, 
and V-R. Color variations should be uncorrelated with the V brightness and with each 
other if the measurement error is 0.03 mag in V instead of 0.015 mag. In either case, the 
Monte Carlo model demonstrates that the fluctuations occur on short timescales. The 
typical observation frequency is ~ 1 day -1 (344 out of 663 observations), 0.5 day -1 (82 
observations), or 0.333 day^ 1 (54 observations). Fluctuations must occur on timescales <^ 1 
day based on the success of the Monte Carlo model in reproducing the light curve. 



3.3. Color Variations 

Figure 3 shows the variations in the optical colors as a function of time. The solid line 
in each of the left-hand panels is the seasonal mean color for the individual data points. 
There is little evidence for a substantial color change in FU Ori during a 0.3 mag decline 
in V. The seasonal means for B-V and V-R are constant to within the photometric errors. 
The variation in the mean U-B color is roughly twice the photometric error but shows no 
obvious trend with time. 

Despite the lack of long-term variation in the colors, there are considerable short-term 
variations. The right hand panels in Figure 4 show the color fluctuations about the seasonal 
means in the left hand panels. The full amplitudes of the color variations are ~ 0.2 mag 
in U-B, and ~ 0.1 mag in B-V and V-R. These variations are 3-5 times the quoted 
photometric errors. 

Figures 4-5 show the correlation of color variations with the V-band fluctuations 
analyzed in §3.2. To construct these plots, we derived separate seasonal means for the 
V-band observations associated with each color observation (because a given color was 
not always obtained with each V measurement) and subtracted the appropriate seasonal 
mean from each data point. The Spearman rank correlation coefficient is 1.3 x 10~ 2 
for 438 S(U-B)SV pairs, 2.6 x 10~ 9 for 626 6(B-V)SV pairs, and 2.9 x 10~ n for 622 
#(V-R)-£V pairs. We fit each correlation with a straight line using the Press et al. (1992) 
subroutine FITEXY, assuming that the la errors in each coordinate are equal to the quoted 
photometric errors. The results yield: 
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5(U -B) = -0.0002 ± 0.0034 - (0.40 ± O.U)SV (5) 

S(B -V) = (-0.24 ± 8.0) x 10~ 4 - (0.12 ± 0.02)5V (6) 

S(V -R) = (-0.015 ± 6.1) x 10~ 4 + (0.15 ± 0.02)SV (7) 

The correlation between <5(U-B) and 5V is weak; U-B becomes redder as the source 
becomes brighter. Both 5 (B-V) and 5 (V-R) correlate well with SV; as the source brightens, 
B-V becomes redder while V-R becomes bluer. 

The slopes of the short-term color variations in equations (5)-(7) differ from the 
apparent slope of the long-period wave in the B and visual light curves. The lack of a clear 
B variation associated with the visual wave suggests S(B — V) ~ C5V, with C*l. This 
behavior suggests that the wave and the short-term fluctuations have different physical 
origins. Analysis of photoelectric data with a longer time baseline is needed to verify this 
point. 

To test the accuracy of the correlation coefficients and measured slopes, we constructed 
artificial color curves using the Monte Carlo model described in §3.2. We matched model 
color curves having amplitudes similar to the observations to model V-band light curves, 
measured the correlation coefficients, and derived the slopes of straight line fits to the 
artificial residuals. We repeated this exercise using known correlations of color index with 
brightness. Random light curves yield no correlation between color index and brightness. 
Correlated changes in color index with brightness yield the measured correlation coefficients 
if the amplitude of the brightness variation is a v = 0.036 ± 0.007, if the slopes and la 
errors of the color variations are those quoted in Equations (5)-(7), and if the photometric 
errors are 0.015 ± 0.01 in V, 0.02 ± 0.01 in B-V and V-R, and 0.07 ± 0.02 in U-B. 

The small color variations make it difficult to find robust correlations among the color 
indices. There is a 3a correlation between U-B and B-V, with a Spearman rank coefficient 
of 4.7 x 10~ 3 . The Spearman rank correlation coefficient is much larger, 0.25, for B-V and 
V-R. To test the importance of these results, we repeated these tests using artificial data 
sets with known correlations between the color indices. The correlation between <5(U-B) 
and 5 (B-V) is always detected in data with small photometric errors; the Spearman rank 
correlation coefficient is 10~ 3 or smaller in all of our trials. The correlation between 5 (B-V) 
and <5 (V-R) has 2a or smaller significance in all of our trials. Reducing the photometric 
errors by factors of 2-3 allows us to recover known correlations at the 3a level. 
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These results provide good evidence that the observed variations in brightness and 
color are intrinsic to FU Ori. The correlations between the B-V or V-R color and the 
V brightness are robust. The measured correlation coefficients are reasonable given the 
magnitude of the color changes and the photometric errors. We next consider the physical 
nature of the brightness changes and then compare these properties with the flickering 
observed in other accreting systems. 



3.4. Physical Nature of Light and Color Variations 

The observed optical light and color variations in FU Ori are small, ~ 0.035 mag in V, 
~ 0.015 mag in U-B, and ~ 0.004 mag in B-V and V-R. To understand the origin of this 
variability, we consider the observed colors as small perturbations about the colors of the 
'average' source in FU Ori. We adopt the mean colors of the system as the average colors, 

[7-5 = 0.84 B-V = 1M V-R=1.15 . (8) 

and correct these colors for interstellar reddening assuming a standard extinction law 
flMathis 1990| ) and Ay = 2.2 mag QKenyon et at. 1988| ), 



(U-B)o = 0.33 (B-V) = 0M (V - R) = 0.60 . (9) 

Figure 6 compares the reddening-corrected color indices of FU Ori with colors for F and G 
supergiants as indicated in the plot. The average colors of FU Ori, shown as the box, are 
offset from standard stellar loci. The solid line in each panel shows how the colors change 
as the source brightens; the colors move away from the supergiant locus if the source fades 
in brightness. Both lines indicate that the best stellar match to the changes in brightness 
and color is a GO supergiant. The dashed lines indicate how the colors change if the slopes 
of the relation between the V-band brightness and the color indices are ±2<r different from 
those derived in Equations (5)-(7). This result shows that the uncertainty in the stellar 
match is less than one spectral subclass. 

The best stellar match to the color variations is correlated with the adopted reddening. 
An F5 supergiant can match the observed color variations for Ay = 3.2 mag; a G5 
supergiant is the best match for Ay = 1.2 mag. Intrinsic colors for Ay ^1.5 mag and 
Ay ^> 3 mag are inconsistent with the G-type optical spectrum. This constraint limits the 
spectral type of the variable source to F7-G3. 
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To find a reasonable explanation for the color variations, we first consider mechanisms 
appropriate for single stars without circumstellar disks. Although we believe an accretion 
disk provides the best model for observations of all FUors, it is important to consider 
alternatives before we identify flickering as the source of the variability in FU Ori. 
Obscuration, pulsation, and rotation are obvious choices for short period variations with 
small amplitudes in an isolated star. Random obscuration events by small, intervening dust 
clouds near the central source are a popular model for rapid variations in Herbig Ae/Be 
stars and other pre-main sequence stars (e.g., |Natta et al. 1997| ; |Rostop china et al. 1997] ; 



and references therein) . The observed pulsational amplitude of the pre-main sequence star 
HR 5999, 0.013 mag ( [Kurtz fc Marang 1995| ), is close to the amplitude observed in FU Ori. 
The period of 5 hr is short compared to the spacing of our data; a similar period in FU Ori 
might well be missed by the analysis described above. Rotational modulations with small 
amplitudes and periods of 1-2 days are observed in many pre-main-sequence stars ( |Kearns 



\et al. 19"9"7| and references therein) and could plausibly produce a similar variation in FU 
Ori. 

Despite its attractiveness, a stellar pulsation in FU Ori seems unlikely. With an 



effective temperature of ~ 6000-6500 K and a luminosity of ~ 200 L (e.g., Kenyon 



1999Q , FU Ori lies within the classical instability strip (see pautschy fc Saio 1995| ). The 



bright Cepheid a UMi (Polaris) has a comparable amplitude and other Cepheids have 
similar periods. Recent period-luminosity relations for Cepheids ( [beast fc Catchpole 1997 



Hindsley fc Bell 1989Q yield periods of ~ 0.8 days for FU Ori if d = 500 pc and A v = 2.2 



mag. Although this period is reasonably consistent with our observations, FU Ori lies well 
above the pre-main-sequence instability strip in the HR diagram ([Marconi fc Palla 1998Q . 



To test the possibility that a pulsational instability might occur anyway, we constructed 
artificial light curves with amplitudes of 0.035 mag and periods of 0.3-2.0 days. We sampled 
these light curves in time as in our real observations and added noise equivalent to our 
photometric uncertainties. Our failure rate for recovering known periods in the artificial 
light curves is small, ^ 10~ 3 , for all periods considered, based on 10,000 artificial light 
curves at each of 20 periods between 0.3 and 2.0 days. The failure rate is largest for periods 
of 0.5, 1.0, and 1.5 days due to the ~ 1 day spacing of the light curve. The failure rate 
decreases to 10~ 4 or less at other periods. Our inability to detect any periodicity in the FU 
Ori light curve and the lack of a theoretical instability strip for the observed luminosity and 
temperature of FU Ori is a strong indication that pulsations do not produce the observed 
variation. 

Stellar rotation is also an unlikely source of the variability. The light curve analysis for 
short periodicities rules out a rotational modulation of the light curve, unless the dark or 
bright spots responsible for the variations vary in size or intensity on timescales of several 
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rotational periods. At the 1-2 day periods that seem most plausible, FU Ori currently 
rotates close to breakup. If the star conserved angular momentum during the rise to 
maximum, the pre-outburst rotational velocity would have exceeded the breakup velocity 
by a factor of ~ 3 (see also Hartmann fc Kenyon 198*5] ) . 



Obscuration events similar to those envisioned in Herbig Ae/Be stars (see [Rostopchina| 
et al. 1997| and references therein) require special circumstances to explain the variability. 



Reddening by small dust clouds requires unusual particles to account for a very steep 
reddening law, Ry ~ 8, that changes sign at V. Light scattered off a reflection nebula 
might account for the different color variations of B-V and V-R, but it is difficult to derive 
as steep a color variation as is observed with simple geometries and a wide range of dust 
properties. Photopolarimetry would test this conclusion. 

Finally, fluctuations in the wind of FU Ori are a less plausible source of brightness 
changes than flickering. Errico et al. (1997) have suggested that the continuum optical 
depth through the wind exceeds unity in the Balmer and Paschen continua. Small variations 
in the optical depth due to inhomogeneities in the outflow might account for small 
brightness changes. If correct, this hypothesis would predict a decrease in the amplitude 
of the variation with decreasing wavelength, because the optical depth in the Paschen 
continuum decreases with decreasing wavelength. For reasonable wind temperatures and 
densities, ~ 5000-10000 K and ~ 10 10 -10 14 cm" 3 flCalvet et al. 1993E |Hartmann fc Calvet 



1995 ), the increase in the optical depth from B to R is roughly a factor of 3 for a gas in 



LTE. There is no evidence for this behavior in FU Ori. 

We conclude that fluctuations associated with a stellar photosphere or wind from the 
disk provide a poor explanation for the rapid, small-amplitude photometric variations in FU 
Ori. The variations in FU Ori have much in common, however, with the flickering observed 
in other accreting systems, such as cataclysmic variables (CVs) and low mass X-ray binary 
systems (Bruch 1992 , 1994; Bruch fc Duschl 1993| ; "Warner 1995| and references therein). In 



most CVs, 0.01-1 mag fluctuations occur on the dynamical timescale, seconds to minutes, 
of the inner disk. Within each flicker, a CV becomes bluer as it gets brighter. The large 
color temperature, ^> 20,000 K, of the variable source also plausibly confines the flickering 
to the inner disk in many CVs. Despite the lack of a good physical model for flickering in 
CVs, it clearly probes physical conditions in the inner disk. 

The observed variations of FU Ori are also plausibly associated with the inner regions 
of a circumstellar accretion disk. The timescale of the variation, <^ 1 day, is close to the 
dynamical timescale of the inner disk, ~ 0.1 day for a 1 M central star with a radius of 4 
R . The temperature of the variable source, ~ 6000 K for a GO supergiant, is comparable 
to the inner disk temperature of 6500 K derived from detailed fits to the spectral energy 
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distribution and the profiles of various absorption lines ( [Kenyon et al. 1988 ; Bell et al. 



1995| ; [lurner et al. 199~7|) . Finally, the amplitude of the variation is similar to that observed 



in other accreting systems. Given this behavior, we believe that the variations in FU Ori 
are flickering and thus provide additional evidence for an accretion disk in this system. 

The main alternative to a variable accretion disk in FU Ori is variations associated with 
a magnetic accretion column. Despite the success of this model in other pre-main- sequence 
stars (e.g., |Gullbring et al. 1996| ), truncating an accretion disk with the large accretion rate, 
~ 1O~ 4 M yr _1 , estimated in FU Ori requires a large magnetic field, ~ 10 kG. Limits on the 
magnetic fields of other pre-main- sequence stars are much smaller, <; 2-3 kG QJohns-Krull 



et al. 1999] ). A much larger field in FU Ori is unlikely. The temperature of the variable 
source in FU Ori, ~ 6000 K, also seems too cool to be associated with a magnetic accretion 
column, where the typical temperature is ~ 10 4 K (e.g., |Lamzin 1998J ; |Oalvet fc Gull bring 



1999). Both of these arguments make a stronger case for flickering as the source of the 



variability in FU Ori. 

To see what we can learn about the inner disk from the variations, we consider several 
simple models for the flickering. We adopt the observed colors - U-B, B-V, and V-R - 
as the colors of the average state of the disk. We assume several sources of variability in 
a steady-state disk, (i) discrete changes in the mass accretion rate, M, through the entire 
disk, (ii) random fluctuations in the flux from any annulus in the disk, and (iii) random 
fluctuations in the flux from specific annuli in the disk. We chose a steady-state disk as 
the average source, because steady disks provide a reasonable fit to the complete spectral 
energy distribution of FU Ori. Several experiments with non-steady disks yield similar 
results. 

Models where the entire disk can vary in brightness do not reproduce the observations. 
We assume a disk composed of discrete annuli with width 5R at a distance R from the 
central star, with SR <C R. Each annulus radiates as a star with the effective temperature 
assigned to the annulus. The arrows in Figure 6 indicate the color variations produced by 
models where the flux from each annulus is a random fraction, between 1.0 and 1.1, of the 
flux from the annulus of a steady-state disk. The color variation of the model clearly fails 
to account for the observed color variation. Allowing all annuli to vary coherently also 
produces color variations that disagree with the observed variation. 

Successful models allow only specific parts of the disk to vary in brightness. If disk 
annuli with the colors of F9-G1 supergiants are the only annuli that vary, the color variation 
of the model follows the slope of the observed variation. 

Our results for flickering in FU Ori are at odds with predictions of the simple 
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steady-state accretion disk model. In steady disk models for FUors, the disk temperature 
rises rapidly from zero at the stellar photosphere, R = R+, to T max ps 6500-7000 K at 
R = 1.36 R* and then decreases radially outward ([Kenyon et al. 1988|). The GO temperature 



of the flickering source has a temperature of ~ 6000 K. Disk material with this temperature 
has R = 2.5 R* in the steady model. Fluctuations in the energy output of this region, and 
the lack of fluctuations in hotter disk material, seem unlikely. If we associate flickering with 
small changes in the mass accretion rate through the disk, M, or in the scale height of the 
disk, H, we expect these to produce larger variations at smaller disk radii. 

Recent calculations indicate that the inner regions of FUor disks may be much different 
than predicted by the simple disk model. The steady-state temperature distribution 
assumes a physically thin disk, H <C R ( [Lynden-Bell fc Pringle 1974| ). FUor disks are 
probably much thicker ( [Lin fc Papaloizou 1985| ; |Ularke et al. 1989] , 1990). Steady models 



that include a self-consistent treatment of the boundary layer between the inner disk 
and the stellar photosphere predict large scale heights, H/R ~ 0.1-0.3 at R ~ 1 — 2 R* 
( Popham et al. 1993J , 1996). Time-dependent models further indicate that H/R can vary 



in a complicated way close to the central star ([lurner et al. 1997 ; Klcy fc Lin 1999|) . Both 



types of model predict that the disk temperature peaks just outside the stellar photosphere 
at R ~ 1.1—1.2 R*. The decline in disk temperature at smaller radii can be as large as 
25%-50%. Applied to FU Ori, these models predict that the disk temperature close to the 
central star is ~ 5000-6000 K, comparable to the temperature derived for the flickering 
source, if the peak temperature at 1.1—1.2 R+ is ~ 7000 K. 

We propose that the flickering source in FU Ori lies between the stellar photosphere 
and the peak temperature in the disk at 1.1-1.2 R+. In the models described above, this 
region produces ~ 5% of the total optical light. Observable variations in the total light 
from the disk - as we have reported here - thus imply significant changes in the physical 
structure of the inner disk. Our data for FU Ori require 50% variations in the light output 
of the inner disk. 

Large changes in the physical structure of the disk can be avoided if the spatially 
thick portion of the disk occults the inner disk. If we view the disk at an inclination i crit 
= ta,n~ 1 H/R, small variations in H/R can produce small changes in brightness and color. 
Rapid variations similar to those observed are possible if H/R varies on the dynamical 
timescale and if the occulted portion of the disk radiates as a GO I star. The required 
change in H/R, ~ 10%, is small compared to the 50% change in light output needed above. 
The required geometry is, however, very special and yields no variation if the real viewing 
angle is much less than i crit . Observations of V1057 Cyg and V1515 Cyg will test this idea, 
because these systems probably have smaller i than FU Ori QHartmann fc Kenyon 1996]) . 
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4. DISCUSSION AND SUMMARY 

Our results provide the first evidence for rapid photometric variations, flickering, 
in a FUor. The amplitude of the variation is small, ~ 0.035, and just detectable with 
photoelectric data covering a long time interval. Observations with smaller photometric 
uncertainties are needed to verify the detection and to place better limits on the color 
variations. Differential photometry using a CCD on a small telescope can achieve the 
required precision, but the field of FU Ori has few bright comparison stars within 15-20 
arcmin. The richness of the field may compensate and allow the high quality photometry 
needed to improve our results. 

Previous attempts to find similar short-term variations in a pre-main- sequence star 
have met with mixed success. Smith et al. (1996) placed upper limits of 0.01 mag on 
short-term fluctuations in four classical T Tauri stars. Gullbring et al. (1996) detected 
flare-like activity in the classical T Tauri star BP Tau, with amplitudes and timescales 
comparable to that observed in FU Ori; Hessman & Guenther (1997) noted similar behavior 
in three classical T Tauri stars, DG Tau, DR Tau, and DI Cep. These studies all interpreted 
the variations with a magnetospheric disk model, where jitter in the magnetically channeled 
flow from the inner disk to the stellar photosphere produces small amplitude 'flares'. We 
prefer to associate the variation in FU Ori with flickering of the inner accretion disk. The 
accretion rate in FU Ori, ~ 10 _4 M Q yr _1 , is a factor of ~ 1000 larger than accretion rates 
derived for T Tauri stars, which makes it difficult to truncate the disk with the modest 
magnetic fields, £ 1-2 kG, detected in pre-main sequence stars. Future observations can 
test the magnetic alternative by placing better limits on any periodic component of the 
photometric variation and by measuring the magnetic field strength. 

However these observational issues are resolved, it is clear that high precision 
photometry can probe the physical conditions of the inner accretion disk of a pre-main- 
sequence star. The amplitudes and timescales of these variations already provide some 
challenge to theory. The amplitude of the FU Ori variation implies large fluctuations in the 
physical structure of the disk on short timescales. Flares and other short-term variations 
in T Tauri stars suggest smaller, but still significant, changes in disk structure close to 
the central star. Recent hydrodynamical calculations show that the disk structure can 
change significantly on longer timescales, but theoretical models do not yet address rapid 
fluctuations in the disk similar to those observed ( [Kley fc Lin 1999] ). Future calculations 



that consider this behavior should lead to a better understanding of mass flow in the inner 
disk in FUors and other types of accreting systems. 
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Table 1 . Optical Photometry 



Julian Date 


V 


B-V 


U-B 


46320.581 


9.28 


1.34 




46325.582 


9.27 






46328.589 


9.27 


1.33 




46357.619 


9.24 


1.32 


0.92 


46378.610 


9.31 






46405.534 


9.32 


1.31 


0.96 


46406.543 


9.33 


1.29 


0.88 


46407.528 


9.28 


1.28 




46466.342 


9.31 


1.32 


0.92 


46703.561 


9.25 


1.36 


1.02 


46710.579 


9.27 


1.34 


0.95 


46761.514 


9.26 


1.30 


0.90 


46773.390 


9.21 


1.34 


0.95 


46774.419 


9.23 


1.33 


0.99 


46831.314 


9.30 


1.32 


0.92 


46852.272 


9.31 


1.33 


1.00 


46855.217 


9.32 


1.33 


1.04 


47183.313 


9.17 


1.36 




47184.308 


9.24 


1.36 


1.02 


47202.361 


9.29 


1.34 


0.96 


47205.233 


9.28 


1.35 


0.99 


47207.319 


9.27 


1.33 




47208.338 


9.29 


1.35 


0.95 


47569.377 


9.32 


1.32 


0.88 


47861.478 


9.42 


1.40 


0.93 


47920.262 


9.39 


1.30 


0.87 


47967.240 


9.38 


1.33 


0.88 


47971.227 


9.42 


1.29 


0.87 


48159.561 


9.45 


1.28 


0.88 


48273.356 


9.42 


1.29 


0.90 


48325.318 


9.46 


1.28 


0.90 


48540.557 


9.45 


1.27 


0.92 


49780.219 


9.44 


1.29 


0.93 


50156.258 


9.45 


1.30 


0.91 


50509.243 


9.54 


1.29 


0.92 


50863.205 


9.57 


1.26 


0.89 
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Table 2. Absorption Indices 



Julian Date CH H/3 Mg I Na I TiOi TiQ 2 
2446722 0.34 0.15 0.07 0.24 -0.01 -0.02 



2447106 


0.23 


0.19 


0.06 


0.23 






2447496 


0.19 


0.16 


0.08 


0.21 






2448404 


0.21 


0.21 


0.14 


0.23 


-0.04 


-0.02 


2449836 


0.33 


0.20 


0.16 


0.26 


-0.03 


-0.06 


2450781 


0.28 


0.22 


0.15 


0.29 


-0.01 


-0.06 


2450983 


0.31 


0.20 


0.19 


0.29 


-0.02 


-0.07 


2451010 


0.31 


0.20 


0.20 


0.29 


-0.02 


-0.09 


2451026 


0.30 


0.24 


0.19 


0.30 


-0.03 


-0.08 
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Table 3. Light Curve Parameters 



Parameter pg + B visual 

mo (mag) 15.55 ± 0.40 

Sm r i se (mag) 5.5 ± 0.4 

T rise (days) 52 ± 25 

to (mag yr" 1 ) 0.0140 ± 0.0015 0.0144 ± 0.0010 
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Fig. 1. — Optical (lower panel) and residual (upper panel) light curves for FU Ori. The blue 
(B) light curve is comprised of 1042 measurements derived from photographic plates and 
photoelectric photometry (Ibragimov 1997 and references therein). The visual light curve 
is ten day means of 7300 visual observations from the AAVSO international database. The 
residuals plot the difference between the B and visual light curves in the lower panel and 
the model light curves described in the text. Both residual light curves show large, short 
timescale fluctuations superimposed on a wave-like variation with a timescale of 3000-3500 
days. 
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Fig. 2. — Photoelectric V-band data for FU Ori. (a) bottom panel: filled circles indicate 
individual data points; the solid line shows seasonal means, (b) top panel: variation of 
individual data points about the seasonal mean light curve. 
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Fig. 3. — Color variations in FU Ori. In each panel on the left, filled circles indicate 
individual data points; the solid line shows seasonal means. Filled circles in the right hand 
panels indicate the variation of each color about the seasonal mean color curves in the left 
hand panels. 
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Fig. 4. — Correlation of U-B variations with V-band variations. The straight line is a 
least squares fit to the plotted points. The error bar in the lower right corner indicates the 
uncertainty of a single observation. 
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Fig. 5. — Correlation of color variations with V-band variations. The lower panel plots the 
variation of B-V with V; the upper panel plots the variation of V-R with V. In each panel, 
the straight line is a least squares fit to the plotted points. The error bars indicate the 
uncertainty of a single observation. 
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Fig. 6. — Flux ratio diagrams for FU Ori. In each panel, a box shows the reddening-corrected 
position of FU Ori, filled circles show positions of F-G supergiants with indicated spectral 
types, the solid line points to the best stellar match to the colors of the variable source, and 
the dashed lines show ±2ct uncertainties in the best stellar match. A GO supergiant provides 
the best match to the observed color variations with an uncertainty of roughly 1 spectral 
subclass. 



